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Abstract. This paper is concerned with a quaslinear parabolic equation including a nonlinear non- 
local initial condition. The problem arises as equilibrium equation in population dynamics with 
nonlinear diffusion. We make use of global bifurcation theory to prove existence of an unbounded 
continuum of positive solutions. 



1. Introduction 

Age-structured equations have a long history (see |21] and the references therein), both for the 
case without and with spatial movement of individuals. In an abstract form, the evolution of an 
age-structured population with quasilinear diffusion is governed by the equations 

d t u +d a u + A(u,a)u + ji(u,a)u = , t>0, aej\{0}, (1.1) 

u(t,0) = / B(u,a)u(t,a)da , t>0 (1.2) 

Jo 

subject to some initial condition at t = 0. Here, the function u = u(t,a) > represents the 
density of a population of individuals with respect to time t > and age a £ /, where / := [0, a m ) 
with maximal age a m . For u and a fixed, A(u,a) is a linear unbounded operator 

A(u,a) : Ei C Eq -> Eo 

defined on some common subspace E\ of an ordered Banach space Eg. It reflects a nonlinear 
spatial dispersion of the individuals. Besides spatial movement, the population undergoes death 
and birth processes with density dependent death modulus ]i(u,a) > and birth modulus 
B(u,a) > 0. Specifically, equation l ll.2|l gives the newborns. 

To understand the asymptotic behavior of the evolution of structured populations, precise in- 
formation about equilibrium (i.e. time-independent) solutions is needed. Clearly, u = is an 
equilibrium solution to (O) , l|1.2|l which therefore has to be singled out in the further analysis. 
Moreover, since u in il.l} , (L2) represents a density, solutions should be nonnegative and thus, 
in our abstract setting, belong to the positive cone Eq of the ordered Banach space Eo- 
For non-diffusive age-structured populations (i.e. A = 0), positive equilibria were established 
e.g. by means of fixed point theorems in conical shells [ 10ll20| . These results were carried over 
to the diffusive case as well [17]. However, the fact that u = is a time-independent solution to 
111. lb . 1 11. 2b allows one also to interpret the problem of finding positive equilibrium solutions as a 
bifurcation problem [7] and so to obtain more insight on the structure of the equilibria set. For 
this we write the birth modulus in the form fi(u,a) = Xb(u,a) and thus introduce a bifurcation 
parameter A > which determines the intensity of the individual's fertility while the qualitative 
structure of the fertility is modeled by the function b. Writing A(«) = A(u, ■) + ]i(u, ■) and 



\£(v)u := / \b(v, a)u(a)da , 
Jo 
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we are concerned in this paper with finding values A > and positive nontrivial solutions 
u : } — > Eq to the nonlinear problem 

d a u + A(u)u=0, aej\{0}, (1.3) 
m(0) = Xl{u)u . (1.4) 

The results derived herein extend our previous results 1161 ITTII on local and global bifurcation. 
In [16 [ it was shown that if A(u) depends sufficiently smooth on u and if for u = 0, the oper- 
ator A(0) possesses maximal Lp-regularity (for a precise definition see the next section), then a 
local curve of positive solutions to ^1.3) . {l.ty bifurcates from the trivial branch (A,u) = (A, 0), 
AeE. This local branch was subsequently extended in fl7[ to a global continuum by apply- 
ing Rabinowitz' global alternative [1 1 J, but for less general diffusion operators. More precisely 
the existence of an unbounded continuum of nontrivial positive solutions (A, u) to jl.3i . jlAl 
was derived under the assumption that the operator A(u) admits a suitable decomposition 
A(u) = Ao + A*(m) with A* being of "lower order". Although the operator A(u) may still 
depend nonlinearly on u, a quasilinear dependence, however, is not covered by the bifurcation 
result of [17J and hence, if A(m) is a second order differential operator, nonlinearities are merely 
allowed in the zero and first order terms. 

The aim of this paper is to remedy this deficiency by establishing a global continuum of positive 
solutions for quasilinear diffusion operators A(m). After recalling (and refining) in Section [2] 
the results of |16J on local bifurcation, we shall show in Section [3] global bifurcation from the 
trivial branch provided some convexity condition (see l|3.2|l for details) holds implying maximal 
Lp-regularity of A(u) for each u. The proof relies on a recent result of Shi & Wang [14J which is 
based on the unilateral global bifurcation techniques of Rabinowitz [11] or rather their interpre- 
tation by Lopez-Gomez |5J. As we shall see then in Section |H the convexity condition d3.2b is not 
necessary provided A and I in {\3\ . (13} depend real analytically on u. Indeed, in this case, the 
analytic bifurcation theory due to Buffoni & Toland [5 ] yields a global smooth curve of positive 
solutions to (|1.4) l. Finally, in Section [5] we give an example of a concrete operator and the 

corresponding function space setting. Further examples to which the results of the present paper 
apply can be found in fl6"l[T7l . Applications of global bifurcation theory to nonlinear population 
models are given e.g. in |18). 

2. Preliminaries 

2.1. General assumptions. If E and F are Banach spaces we write £(E, F) for the set of linear 
bounded operators from E to F, and we put C{E) := C(E,E). The subset thereof consisting 
of compact operators is denoted by /C(E, F) and /C(E), respectively. Isom (E,F) stands for the 
set of topological isomorphisms E — > F. By E "— > F we mean that E is compactly embedded in F. 

Throughout the paper we assume that Eq is a real Banach space ordered by a closed convex cone 
Eq" and Ei is an embedded Banach space such that 

the embedding E\ » Eq is dense and compact . (2-1) 

We fix p 6 (1, oo ) and set E e := (Eq, Ei) s ,p with (•, being the real interpolation functor 
for g := g(p) := 1 — 1/p. For each 8 6 (0,1) \ {1 — 1/p} we let (■, -) e denote an admissible 
interpolation functor, that is, an interpolation functor (■, - )e such that the embedding 

Ei Ee : = (Eo/Ei)e 

is dense. Note that the embedding Eg "^-> E# is compact for < # < 9 < 1 (see [3, 1.Thm.2.11.1]). 
The interpolation spaces Eg, < 9 < 1 are given their natural order induced by the cone 
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:= E e n Eq". We suppose that 

int(E+)^0. (2.2) 

Recall that a m £ (0,oo] and set / := [0,a m ). Observe that a m = oo is explicitly allowed. We 
introduce the spaces 

E := L P (/,E ) , Ex := L p (/,Ei) D Wj(/,E ) 
and recall the embedding 

Ei ^BUC(J,E g ), (2.3) 
where BUC stands for the bounded and uniformly continuous functions. Thus, the trace 

y u := u(Q) , mGEj, 

yields a well-defined operator 70 £ £{E lr E g ). We let E+ := L+(/,Ei) n W^(J,E ) denote the 
positive cone of Ej and put E^ := E+ \ {0}. We fix a Banach space F such that 

Ei F (2.4) 

and let E denote an open connected zero-neighborhood in F. Then Ei := En Ej is an open 
connected zero-neighborhood in E^. Suppose that for some d £ we hav43 

A £ C 1 (E,£(Ei,E )) , <eC 1 (E,£(Ei ( £»)) (2.5) 

and that for each m £ Ei, the operator A(u) possesses maximal L p -regularity, that is, 

(d a + A(«),7o) £ Isom(E 1 ,E x E c ) , u £ Ej . (2.6) 

Then 

(u ^ T[u] := (3 fl + A(m), 7o )" 1 ) £ C(Ei,£(E x E c ,E a )) (2.7) 
due to continuity of the inversion map B 1— > B _1 . We suppose that 

T[«](CV) e£+(E C/ Ei), HGEi, (2.8) 
that is, T[u] (0, ■) maps the positive cone E+ into the positive cone Ej~. Set 

Q(u) :=*(k)T[u](0,-), «e Ei, 
and note that from j2.5\ and E# » E e we have 

Q £ C(Ei,£(E c ,E tf )) nC(E 1 ,/C(E G )) . (2.9) 

We further suppose that 

Aq" 1 > is a simple eigenvalue of Q(0) with eigenvector <I>o G m t(E^), 

and there is no other eigenvalue of Q(0) with an eigenvector in E^ . 

For our global bifurcation results we shall strengthen some of the conditions. We point out that 
assumptions above are met in many applications, roughly i2.7i is satisfied e.g. for second order 
operators in divergence form and | |2.8t is the maximum principle for parabolic equations. We 
shall be more specific in Section [5] Regarding assumptions ||2.4| | and | |2.10| | we note: 

Remarks 2.1. (a) Let a m < 00. If a £ [0, 1) and s £ [0, 1 - a), then Ei W*(J, E a ). 

Proof. This follows from a generalized Aubin-Dubinskii lemma, see [4[ Thm.1.1]. □ 

(b) 1/ Q(0) £ /C(£ g ) is strongly positive, i.e. if Q(0)O £ int(E+)/or eac/z O £ E+ \ {0}, then (f27TOT > 
/toWs wif/j A^" 1 equals the spectral radius o/Q(0). 

Proof. This is a consequence of the Krein-Rutman theorem [8. Thm. 12.4]. □ 



^Observe that this notation includes that A = A(u,a) in jl.lt depends in a local way on age a. 
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Let us note that (A, m) e K x solves 11.3b , 11.4b if and only if 

K = T[tt](0,tt(0)), tt(0) = AQ(u)u(0) , (2.11) 

which follows by plugging the solution of 1 11. 3b , given by the first identity of 12.11b , into 11.4b . 
Equivalently setting S := T[0] we see that (A, u) G JR. x Zi solves 11.31 , 11.4b if and only if 

h = S((A(0) - A(u))u, X£(u)u) . 

We shall use both characterizations of solutions in the sequel. Note that 12. lit implies that u(0) 
(if nonzero) is an eigenvector of Q(u) with eigenvalue A -1 . 

2.2. Local Bifurcation. By what we have just observed, solving 11.31 , 11.41 is equivalent to finding 
the zeros (A, u) of the function F : R x Zj — > Ej defined as 

F(A,m) := u - S((A(0) - A(u))u, M(u)u) , (A,«) G R x Z t . 

Let 

6 := {(A, u) G IR x Z x ; F(A, «) = 0} . 
Clearly, (A, u) = (A, 0) for A 6 R gives a trivial branch in & of solutions to 11.31 1, 11.4b . We shall 
next show that a nontrivial branch of positive solutions bifurcates from this branch at the point 
(A, u) = (Ao, 0). For this we first show that the Frechet derivative F„(A, u) with respect to u, 
given by 

F u (X,u)[(p} =0-S((A(O)-A(u))0,Ai(u)0) - S( - A u (u)[<p]u, A.£ u (u)[<p]u) , <p G E x , 

(2.12) 

is an index zero Fredholm operator. This will follow from the following observation: 
Proposition 2.2. Lef (A, u) G R x Zj be fixed and set 

T<p:= .F(A,m)4>:=4>-S((A(0)- A(m))</>, A£(u)<p) , <p e E a . 
T/xew T G £(Ei) is a Fredholm operator of index zero. More precisely, 

ker(7") = {F[ U ](0,ro);n;Gker(l-AQ( M ))} 

and 

rg(^) = {heE 1 ;h{0)+M(u)T[u](d a h + A(0)h,0) G rg(l - AQ(m)) } (2.13) 

dim(ker(J r )) = codim(rg(J r )) = dim(ker(l — AQ(m))) < oo . 

Proof. The idea of the proof is the same as in (16l Lem.2.1] and it is rather the functional analytic 
setting that has to be modified slightly. For the reader's ease we include the proof here: By 
definition of S = T[0], the equation J-<p = h for <p, h G E^ is equivalent to 

d a( p + A(u)<p = d a h + A(0)h, (2.14) 

<p(0) - U(u)<p = h(0) . (2.15) 

From 12.14|l it follows 

cp = T[u] (d a h + A(0)«, 0) + T[u] (0, (p(0)) (2.16) 
and when plugged into 12.15b we obtain 

(l-AQ(u))cp(0) = h(0)+M(u)T[u}(dah + A(0)h,0) . (2.17) 

The statement of the proposition is trivial if A = 0, so let A ^ 0. If 1 / A belongs to the resolvent 
set of Q(u) G fC(E g ), then d2TT6b , j27T7b entail a trivial kernel ker(7") = {0}. Moreover, in this 
case, for an arbitrary h G Ej, there is a unique <^>(0) G Eg solving 12.17b as its right hand side 
belongs to E 5 . Consequently the corresponding <p G Ej given by 12.16b is the unique solution to 
J-<p = h. This gives the assertion in this case. 
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Otherwise, if 1/A is an eigenvalue of Q(u) E /C(E S ), then | |2.16| |, j2.17|l yield the characteriza- 
tion of ker(J r ) and rg(J-~) as claimed. In particular, since T[u) is an isomorphism, we deduce 
dim(ker(J r )) = dim(ker(l — AQ(w))) which is a finite number because 1/A is an eigenvalue of 
the compact operator Q(u). Moreover, rg(J r ) is closed in Ei since M := rg(l — AQ(w)) is closed 
by the compactness of AQ(m) and due to i2.3h i2.5i , and (2.7> . Next, to compute codim(rg(J 7 )) 
note that 

codim(M) = dim(ker(l - AQ(m))) < oo , 

hence M is complemented in E £ which yields a direct sum decomposition Eg = M © N. Denoting 
by Pm £ jC(E e ) a projection onto M along N, we set 

P/z := S(d a h + A(0)h, P M h(0) - (1 - P M )M(u)T[u] (d a h + A(0)h,0)) , keE 1( (2.18) 

and obtain P 6 £(Ei) from O, (23), and (27). Since 

(3. + A(0)) (P/i) = 3 fl ft + A(0)/z , 7o(P^) = *Wz(0) - (1 - P m )A€(m)T[m] (3«ft + A(0)h, 0) , 

the characterization 112.13b actually implies that P maps Ej into rg(J-~). Furthermore, if h E rg(J r ), 
then | |2.13l l also ensures 

¥h = S(d a h + A(0)h,h(0)) = h, 
so P(rg(J r )) = rg(J r ). Thus P 2 = P with rg(P) = rg(J r ) is a projection and 



Ei = ker(P) © rg(J r ) . (2.19) 
Since S is an isomorphism, we obtain 

ker(P) = {lieE i; d a h + A(0)h = 0, fe(0) 6 N} , (2.20) 

from which we deduce dim(ker(P)) = dim(N) and the statement follows. □ 



For future purposes let us explicitly state the following decomposition of E^. 
Remark 2.3. The direct sum decomposition 

Ex = span(S(0,cl> )) © rg(F M (A ,0)) 

holds. 

Proof. Taking (A, u) = (Ao,0) in Proposition 12.21 and noticing that T(Xq, 0) = P«(Ao,0), the 
assertion follows from (219)1 and (270) with N = ker(l - A Q(0)) = R • 4> . □ 

Corollary 2.4. For each (A, u) E R x I4, f/ie Frechet derivative F u (\,u) E £(Ei) is a Fredholm 
operator of index zero. 

Proof. Set 

K{u)<p:=S(-A u {u)l<p]u,M u {u)l<p}u) , <p E E x . 

Then, by (2.41 and (2.51 1, K(u) coincides with the Frechet derivative 

K(u) = D W S( - A{w)u, M(w)u)\ w=u E £(F,Ej) C /C(E X ) . 

Consequently, Proposition 12.21 and (2.121 show that F U (A, u) is a compact perturbation of a Fred- 
holm operator of index zero, so F U (A, u) itself is a Fredholm operator of index zero. □ 

Next, we verify that we may apply the Crandall-Rabinowitz theorem on local bifurcation for the 
map F. 

Corollary 2.5. The kernel of F u (Aq,0) is one-dimensional, i.e. ker(F u (Ao,0)) = R ■ S(0,<&o)> an d the 
transversality condition F\ /U (Aq, 0) [l, S(0, <3>o)l r g(Eu(Ao, 0)) *' s satisfied. 
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Proof. It readily follows from j27TOl >, $2A2\ , and Proposition|Z2]that ker(F„(A 0/ 0)) = R • S(0,O ). 
Moreover, (|2l0t and EU l imply F A/M (A ,0) [l,S(0,4> )] = -(0, A Q ^q). Suppose that 

-(0,A - x «>o) Grg(F M (A ,0)) . 
Then Aq^Oq G rg(l — AqQ(O)) due to Proposition 12. 121 contradicting the fact that 

rg(l-A Q(0))nker(l-A Q(0)) = {0} 
since A^ 1 is a simple eigenvalue of Q(0) according to | |2.10|| . □ 

Based on the foregoing observations we are in a position to apply the celebrated Crandall- 
Rabinowitz theorem [6] on local bifurcation and obtain a branch in & of positive solutions to 
CQK CO)- The following result has been observed in [16 1. 

Theorem 2.6. Suppose |(2J]|, |(Z2|> , (ZM . O, $ZM> and <EH0j) hold. Then there are e > and a 
continuous function (A, «) : (— e, e) ->RxSi smc/z t/zaf the curves 

■ = {(X(t),u(t)) ;0<±t<e}ce 

bifurcate from the trivial branch {(A,0) ; A G R} af (A(0),u(0)) = (Ao,0) and 

w(f) = fS(0,O ) + o(f) ast^O. (2.21) 

Near the bifurcation point (Aq, 0), a/Z nontrivial zeros ofF lie on the curve U & + . Moreover, 

•S + \{(A ,0)} C (0, M )xt{ 

and ST n (0,oo) x E+ = 0. 

Proof. According to Corollary 12.41 Corollary 12.51 and [6[ we only have to prove the statements on 
positivity of the last sentence. From (|2.21l l and l|2.10b it follows that 

f- 1 70 M(f) = <J> + 7o^p e int(E+) as t -> , 
whence 7ofi(f) 6 E+ provided t G (0,e) is sufficiently small. Since 

a(0 = r[«(0](o,7o*(0) / 

we conclude il(t) G from l!2.81 l and A(f) > for t G (0, e) from the assumption Ao > 0. 
Consequently, & + \ {(A ,0)} C (0, oo) x E^ . The same argument shows that Mr D (0, oo) x E^ 
is empty. □ 

Note that if A and i in 112.61 1 are real analytic, then so is the curve (A, u) : (— e, e) — > 1R x Li since 
F : 1R x Li — > Ei is real analytic in this case. 

The direction of bifurcation may be determined from the second relation J2.11I I. For instance, if 
Q(w) G K.(Er) is strongly positive for u G Lj, then Il2.11b implies Ar(Q(w)) = 1 for any positive 
solution (A, u) G IR x Ej of ( 11.3b , d 1 ,4b according to the Krein-Rutman theorem. Conditions may 
then be imposed which ensure r(Q(u)) < 1 so that bifurcation is necessarily supercritical. For 
further details we refer to 11161 . 

We next prove a compactness property of the solution set & that we shall use in the coming 
subsections. 

Lemma 2.7. Any bounded and closed subset of& is compact in ExEj. 
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Proof. Let (A n ,M n ) ne ]N be any sequence in a closed subset of & such that || (A n , Mn) HkxEi < c 
for all n E N and some cq 6 1R. Then 

u n = T[m„](0,u„(0)) , k„(0) = A„Q(m„)m„(0) . (2.22) 

Due to l|2.4|l , we may assume without loss of generality that (A n , u n ) — > (A, m) in ]R x E. Ac- 
cording to and this implies 7>„](0,-) -> T[w](0,-) in £(E 6 ,Ei) and Q(m„) -> Q(u) 
in £(E e , E#). Also note that | |2.3l l entails ||mh(0)||e < c for neN. Consequently, from j2.22i we 
derive 

IK(0)||e, < |A„| \\Q(un)\\c(E g ,E g ) II"«(0)||e 6 < c , n e N , (2.23) 

and we thus may assume without loss of generality that M n (0) — > in E e since E# "->• E g . 
Now shows that » = AQ(m)u and u„ = T[%](0,%(0)) -> r[w](0,t7) in E\. Clearly, 

u = T[u] (0, v) since u n — > m in F and we conclude that (A n , w n ) — > (A, u) in ]R x Ej. This proves 
the assertion. □ 



3. Global Continua 

We shall show that the local curve ^ + provided by Theorem 12.61 is contained in a global con- 
tinuum of positive solutions to il.3\ , HAl . In this section we make use of the unilateral global 
bifurcation theory in the spirit of Rabinowitz's alternative IITT1I91 as proposed in [14J. In Section|4] 
we shall give a slightly different approach by means of analytic bifurcation theory [ 5 1 . 
In order to apply the results of [14J we have to strengthen certain conditions. More precisely, in 
the following we suppose in addition to the assumptions stated in Section [2] that 

Eq and Ei are separable (3-1) 

and we strengthen assumption j2.6\ to a convexity condition 

(d a + [(1 -a)A(O) +oA(u)],7o) G Isom(E 1 ,E x E g ) , hGEj, a e [0,1] . (3.2) 

We also assume that 

i(0) e £+(E lr E e ) (3.3) 

and 

Q(0) e K.(E e ) is strongly positive . (3.4) 

Recall that this last assumption implies | |2.10| | with A^" 1 being given by the spectral radius of 
Q(0) according to Remark 12.11 Note that Eq is separable since E\ is separable and dense in Eq. 
Also note that if £o is reflexive (and separable), then Eg is separable. 

Assumption (|3.1|l implies the following remark which is needed in order to apply the result 
of fl4l 

Lemma 3.1. Ej can be equipped with an equivalent norm which is differentiable at any point different 
from 0. 

Proof. Due to fl3l . as Ei is separable since E\ is, the statement that Ei can be equipped with 
an equivalent norm which is differentiable at any point different from is equivalent to saying 
that the dual space E^ is separable. But since E' is separable, the space E Q = L„i(J,E' ) (with 
1 = 1/p + 1/p') is also separable and, moreover, densely injected in E^ since E^ is densely 
injected in Eq. Whence E^ is separable. □ 



The convexity condition l|3.2l l yields: 
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Lemma 3.2. For (A, «)eRx and a G [0, 1], £/ze operator 

(l-a)F„(A,0) + aF M (A,«) e £(E X ) 

is Fredholm of index zero. 



Proof. Let (A,m) elxEj and a e [0,1]. Noticing that 

(l-a)f„(A,O)[0]+aF M (A,u)[4>] = (p - S(«(A(0) - A{u))<p, X[(l - ot)l{0) + od{u)]<p) 

-dlS{- A u (u) [<p] u , U u (u) [<p] u) 
for <p G Ei, the proof is the same as in Proposition 12.21 and Corollary 12.41 by taking | |3.2| | into 



Now we can prove that there is global continuum of positive solutions to (|1.3|l . jlAl . 

Theorem 3.3. Suppose <|2J}, ||Z2), (fZUl, (23), (23)l . and <J3TT])— d3T4j ftoZd. T/ien ffrere z's a connected 
component € + of& containing the branch & + such that € + \ {(Ao, 0)} C (0, oo) x Bit. Moreover, one 
of the alternatives 

(i) € + is unbounded in E x Ei, or 

(ii) £ + intersects with the boundary E x 3Ei 

occurs. In particular, ifL — W, then € + is unbounded in RxEj. 

Proof. Due to Lemma 13.11 and Lemma 13.21 we may apply fl4l Thm.4.4] and deduce that M. + is 
contained in a connected component £ + of & and one of the alternatives 

(a) £ + is not compact in E x Ei, or 

(b) <t + contains a point (A*, 0) with A* 7^ Ag, or 

(c) £ + contains a point (A, z) with z^0 and z 6 rg(F„(Ao,0)) 

occurs, where we have used in (c) the decomposition Ei = span(S(0, <J>o)) © r g(fw(^0/0)) stated 
in Remark 12.31 Notice that, owing to Lemma [2.71 (see [14, Rem.4.2]), alternative (a) is equivalent 
to saying that 

(i) £ + is unbounded in E x Ei, or 

(ii) £ + intersects with the boundary E x d'Li . 

According to Theorem l2.61 the component € + near the bifurcation point (Ao, 0) coincides with 
Next, we show that <£ + \ {(Ao,0)} C (0, 00) x E+. Indeed, if £+ leaves (0, 00) x E^ at some 
point (A,w) 6 £ + fl R x Ei with (A,w) i (0,oo) x E+, there is a sequence ((Ay,My)); € N in 
£+ fl (0,oo) x E+ such that (A ; -,w ; -) ->• (A,m) in E x E a . Clearly, A > and m e E+ with A = 
or m = 0. But since (A, u) £ 6, we readily deduce from (J2JJ} that A = implies u = 0. Hence 
« = in any case, i.e. (Ay, uf) — > (A,0) in R x E^. Again by 112.111 1, we have 



of (vj) (which we do not index) which converges to some v in F. From Il2.7b and (|2.9l l we deduce 
T[uj] T[0) = Sin £(E x E^Ej) and Q(«y) Q(0) in £(E g ,E#). As in E23t we then obtain 
from and (23J that 



from which we conclude that the sequence (vj(0)) is bounded in E$ Eg. So, extracting a 



account. 



□ 




W;'(0)||e„ < c||w/(0)|| E{ < c||My|| El , 




E+ . Letting j — > 00 
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from which we first deduce that A > since otherwise w = implying the contradiction v = 0. 
Consequently, w £ is an eigenvector of Q(0) to the eigenvalue 1/A. Thus A = Aq and 
zo = <x.<3>o for some a > according to l|2.10l l, hence (A, u) = (Ao, 0). Therefore, € + leaves the 
set (0, oo ) x E+ only at (Ao,0) and £ + \ {(Ao, 0)} is thus contained in (0, oo) x E^. In particular, 
alternative (b) above does not occur. We finally show that alternative (c) does not occur as well. 
Suppose to the contrary that £ + contains a point (A, z) with z ^ and z = F u (Xq, 0)£ for some 
I £ Ei. Then z £ E+ and £ - z = S(0, A ^(0)O- Recall from Corollary ES that 0* := S(0,O ) 
with <&o from d2.10b satisfies <p* = S(0, A £(0)^*)- Since O 6 int(E+), we find k > such 
that kO + £(0) -z(0) belongs to E+. Set t/> := k<£* + £ - z. Due to ip = \ S(0,£{0){ip + z)) we 
conclude 9 fl i|> + A.(0)ip = on the one hand from which tp = S (0, t/?(0)), and 

xp(0) = Ao£(0)i/> + A £(0)z 

on the other. Combining these two observations we derive the equation 

(l-A Q(0))!/>(0) = A £(0)z. (3.6) 

However, since A,^ 1 is the spectral radius of the strongly positive compact operator Q(0) and 
since Ag£(0)z £ £+ by (|3.4l l, equation | |3.6| | has no positive solution according to |8, Cor.12.4] 
contradicting ip(Q) = k$o + £(0) — z (0) £ • Therefore, alternative (c) above is impossible and 
the theorem is proven. □ 

4. Global Branches in the Analytic Case 

In this section we shall show that assumptions I |3.1| |-| |3.4| | of Theorem |3.3| are not needed to extend 
the function (A, u) from (0, e) to (0, oo) provided that A and I are real analytic. In this case we 
obtain a slightly better result. So, let us strengthen assumption l|2.5|l to 

A £ C w (2:,£(Ei,E )) , £ e C a; (E,£(Ei,E I? )) . (4.1) 

As noted in Section[2j the function (A, u) is real analytic in this case. 

Theorem 4.1. Suppose lEJ), lO, lO, J231> . fZM , < f27TOl >, and igj] froZd. T/ien f/zere is a contin- 
uous curve 9i + = {(X(t),u(t)) ; f £ [0,oo)} c 6 extending & + . The curve d\ + \ {(Ag,0)} lies in 
(0, oo) x E^ and has at each point a local analytic and injective reparametrization. Moreover, one of the 
alternatives 

(i) || (A(f), a(t) ) || RxEl oo as t or 

(ii) u(t) — > 3Si as t ^ oo 

occurs. In particular, if E = F, f/zen (f) occurs. 

Proof. Due to Corollary 12.41 Corollary 12.51 and Lemma |2~71 we may apply [51 Thm.9.1.1]. Con- 
sequently, there is a continuous curve £H + = {(A(f),w(f)) ; t £ [0, oo)} c 6 extending & + and 
having a local analytic and injective reparametrization. For this curve $H + , one of the alternatives 

(i) || (A(f), u{t)) HrxEj -> oo as £ or 

(ii) u{t) — > 3Si as f — >■ oo, or 

(iii) 5H + is a closed loop, i.e. there is a minimal T > such that D\ + = {(A(£),z7(f)) ; < t < t} 
and (A(t),«(t)) = (A(0),u(0)) = (A ,0), 

occurs. Moreover, if (A(£i), fi(ii)) = (A(f2), "(£2)) ror some £1 7^ £2 with 

ker(F H (A(h),M(fi))) = {0}, 
then (iii) occurs and t\ — £2 | is an integer multiple of t. Finally, the set 

{f >0;ker(F M (A(f),w(f))) ^{0}} 
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has no accumulation points. So, the assertion follows provided we can prove that the curve 
$K + \ {(Ao,0)} lies in (0, oo) x Ej~ and that alternative (iii) does not occur. For this we use an 
argument similar to [5, Thm.9.2.2]: Set 

T* := sup{T > 0; (A(f),fi(f)) £ (0,oo) x E+ for < t < T} 
and note that T* > £ according to Theorem 12.61 Assuming T* < oo, there is a sequence (fy) with 
tj T* such that (Aj,Uj) := (A(f,-),fl(fy)) £ (0,oo) x E^ converges to (A,m) := (A(T*),fl(T*)) in 
R x Ei and (A, w) ^ (0, oo) x E^. But then, the same argument as in the proof of Theorem 13.31 
yields (A(T*),w(T*)) = (Ao,0). Since Theorem 12.61 implies that the bifurcation curve which lies 
in R + x E^ and passes through (Aq, 0) is near this point uniquely determined by we derive 
that (A(£), «(f)) belongs to M, + for f less, but close to T*. Consequently, there are sequences r k \ 
and sjt \ such that (X(r k ),u(r k )) = (A(T* - s k ),u(T* - s k )) and ker(F u (X(r k ),u(r k ))) = {0}. 
But then, as stated above, the minimally chosen t > with (Ao, 0) = (A(t),w(t)) divides 
— s k — r k for each k £ N which is obviously impossible. Therefore, = oo and alternative 
(iii) above does not occur. □ 



Note that it is not claimed in Theorem 14.1 1 that VV~ is a maximal connected subset of ©. Other 
curves or manifolds in © may intersect We also point out that alternative (i) in Theorem 14. II 
is stronger than saying that 9\ + is unbounded in R x Ei (see Theorem I3.3I I. 

5. Example 

Let us consider an example. Let u = u{a,x) denote the distribution density of individuals of 
a population with age a £ / := [0, a m ) at spatial position x in a bounded, sufficiently smooth 
space region Q C R", where a m £ (0, oo) denotes the maximal age. Suppose that the individuals 
move within Q and that dispersal speed d > depends on the local overall population; that is, 
suppose that movement is governed by a density-dependent diffusion term — div x (d(U) V x m), 
where 

U(x) := J u(a,x)da (5.1) 

is the overall population at spatial position x £ Cl. Assume further that individuals cannot 
leave the space region O so that the behavior on the boundary dCl is described by a Neumann 
condition d v u = 0, with v denoting the outward unit normal to dCl. For a given distribution u let 
ji{U, ■) and Ab(U, •) denote the death rate respectively the parameter-dependent birth rate. Then 
equilibrium (i.e., time-independent) solutions to the corresponding evolution problem satisfy the 
following equations: 

d a u - div x (rf(Lf(x))V x w) + }i(U{x),a)u = 0, a £ (0,a m ) , xed, (5.2) 

ram 

u(0,x) — A/ b(U(x),a) u(a,x) da, x£Q, (5.3) 

Jo 

d v u(a,x)=0, ae(0,a m ), x£8Q. (5.4) 

Fix p £ (ft + 2, oo) and set 

Ei := W£ N (n) := {v £ W^(Q);d v v = 0} «-+ E := L p (Q) 

and 

Ei := L p (/,W2 N (n)) n Wi(7,L p (0)) , E := L p (/,L p (Q)) . 

Clearly, l|2.1b and | |3.1| l hold. Observe then that the interpolation result [TBI Thm. 4.3.3] and 
Sobolev's embedding theorem imply 

E 1 = W p 2 N (H) ^ E e := (L v {n),Wl N ) x _ yvv = W^ 1/p) C\tt) . (5.5) 
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Thus int(E+) / while Remark |2J1 implies that 

E 1 ~+V:=W p (J,E i ,) 

for some d > 1 - 1/p = g and some s £ (0,1/p), where E e = W^(Q). Let d £ C 4 (R) satisfy 

d(z) > d > Oforz £ Rand let ji,b eC 4 (Rx [0,%J,JR+) withfc(0,-) / 0. For u 6 F ^ Li (/,E#), 
set 

/•«m 

U := u(a, Ada £ £# 

and define 

A(u,a)w := — div x (d(!J) V^z^) + ]i(U,a)w , w £ Ei , u £ F , aEj. 
Then, by OH Prop.4.1], 

a e c 1 (F ; L 00 (/ / £(w2 N (n) / L p (n)))) 

and thus in particular A G C 1 (F,£(Ei,Eo)). Since l|T9l Prop .4. 11 together with the multiplication 
result of Thm.4.1] ensure leC 1 (F, £(Ei, £#)) as well with 

M(v)u:— / Ab(w,fl)u(fl)dfl , c £ F , ksEj, 
jo 

condition (|2.5b holds. Note that for a £ [0, 1], w £ F and a»eEi we have 
A a (u, -)w : = (1 - a)A(0, -)w + uA(u, -)w 

= -div x ([{l - a)d(0) + ocd(U)]V x w) + [(1 - a)pi{0, •) + xji(U, -)}w 

with (1 — a)d(0) + ad(U) > d. Hence, for a £ [0,1], u £ F, and a £ / the operator —A a (u,a) 
is resolvent positive, generates a contraction semigroup on each Lq(Q), 1 < q < 00 (see U), 
and is self-adjoint in L 2 (C1). Hence III.Ex.4.7.3,III.Thm.4.10.10] entail (32). Since for w £ F 
fixed, the mapping A(u, ■) : [0,fl m ] — > £(Wp N (Q),Lp(Q)) is Holder continuous, there is a 
unique positive evolution operator Yl u (a,<j), < cr < a < a m on £0 corresponding to A(«, •), 
see |3l II.Cor.4.4.2.,II.Thm.6.4.2]. In particular, T[u](0, •) = n„(-,0) G £+(E c ,Ei) for w £ F and 
l{0) G £+(Ei,E G ), that is, (EH) and j33) hold. Also note that the maximum principle ensures 
that rio(fl,0) G /C(E C ) is strongly positive for each a e J \ {0}. As b(0, ■) / Owe conclude (see 
CH Sect.3]) that 

/•Am 

Q(0) = ] b(0,fl)n (a,0)dfle/C(E g ) 

is strongly positive, whence (|3.4b . Consequently, we are in a position to apply Theorem 13.31 
and deduce that there is an unbounded continuum of positive solutions (A, u) in (0, 00) x E^ 

to E3-d5p. 

Let us point out that this is just one simple example which can be extended in various ways. For 
instance, one may consider more general (uniformly elliptic) differential operators that depend 
also locally on age a. Also the regularity assumptions are not chosen optimally and the phase 
space Eq can be any Lq(Cl) provided q £ (1, 00). In Theorem 14 . 1 1 one may also take Eq = Li(Q), 
where one may check the analyticity condition (|4.1) l with the help of [12, 5.Thm.4]. We refrain 
from giving details and refer to llT6l H7 | for other examples. For more concrete applications of 
global bifurcation results we refer e.g. to [18J and to forthcoming research. 
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